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Abstract 

This paper concerns time-dependent scattering theory and in particular the concept of time 
delay for a class of one-dimensional anisotropic quantum systems. These systems are described by 
a Schrodinger Hamiltonian H = —A + V with a potential V{x) converging to different limits Vi 
and V r as x — > — oo and x — > +oo respectively. Due to the anisotropy they exhibit a two-channel 
structure. We first establish the existence and properties of the channel wave and scattering 
operators by using the modern Mourre approach. We then use scattering theory to show the 
identity of two apparently different representations of time delay. The first one is defined in terms 
of sojourn times while the second one is given by the Eisenbud-Wigner operator. The identity 
of these representations is well known for systems where V(x) vanishes as |a;| — > oo (Vi = V r ). 
We show that it remains true in the anisotropic case Vi ^ V r , i.e., we prove the existence of the 
time-dependent representation of time delay and its equality with the time-independent Eisenbud- 
Wigner representation. Finally we use this identity to give a time-dependent interpretation of the 
Eisenbud-Wigner expression which is commonly used for time delay in the literature. 
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I. INTRODUCTION 



Time delay is an important concept in scattering theory. In the simplest situations 
it expresses the excess time that scattered particles spend in the scattering region when 
compared to free particles subject to the same initial conditions. A positive time delay 
means that particles take more time to pass through the region where they are influenced 
by the interaction than particles propagating freely through the same region. A negative 
time delay means that on average the scattered particles are accelerated by the effects of 
the interaction. 

Different approaches to the definition of time delay and related concepts, together with 
various applications to physical problems and a considerable number of references, have been 
presented in a recent review on time delay by de Carvalho and Nussenzveig [lj . In quantum- 
mechanical scattering theory time delay is relevant in particular for the characterization 
of resonances ( 3], [3], ^J), and it enters Levinson's theorem relating scattering data to the 
number of bound states (111, 191). Time delay is related to the density of states in mesoscopic 

n n n q 

conductors p), to the virial coefficients in statistical mechanics p|), and it plays 
a role in the study of chaos ([jj, Q]). 

It is clear that the natural framework for defining the notion of time delay is that of time- 
dependent scattering theory. We refer to the book l| for a general account of this theory 
and just recall its basic ideas. One considers a physical system described by a Hamiltonian 
H acting in some Hilbert space 7i. The following questions are then considered: (i) Given 
a state vector ip G 7i at time t — 0, what kind of asymptotic behavior can one expect to 
see for ip t = e~ lHt ip as t — »• ±oo ? Often it is possible to show the absence of singular 
continuous spectrum of H, meaning roughly that - except for admitting superpositions - 
only two types of asymptotic behavior are possible: ip t may stay essentially localized in a 
bounded region of configuration space (bound state) or disappear from each bounded region 
(scattering state). For state vectors of the second type one then asks: (ii) Can their evolution 
be described asymptotically (i.e., for large values of \t\) in terms of a simpler Hamiltonian 
Hq called the free Hamiltonian (or more generally in terms of a family of free Hamiltonians 
called channel Hamiltonians) ? The answer is positive in particular in situations where one 
can verify the existence of the M0ller wave operators; this then leads to the introduction 
of the scattering operator S that establishes the link between incoming and outgoing free 
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states. A final important point then is: (iii) Can one describe the evolution of all scattering 
states in terms of H (or of the family of channel Hamiltonians) ? If so the theory is said to 
be asymptotically complete. Asymptotic completeness implies in particular the unitarity of 
S. Once the preceding dynamical questions have been settled one may introduce physical 
quantities like time delay or scattering cross sections and determine their properties and 
their relation to the S'-operator. 

Very detailed results are known in time-dependent scattering theory for Hamiltonians 
describing a single non-relativistic particle in a potential V(x) in n-dimensional space under 
the assumption that the potential tends to zero at large distances (i.e., as |x| — > oo). On the 
other hand the literature on scattering theory in highly anisotropic situations, for example 
with potentials assuming different limits in different directions, is rather sparse (except for 
scattering relative to a periodic Hamiltonian and for one-dimensional Hamiltonians to be 
discussed below). For n > 2 the points (i)-(iii) have recently been investigated for potentials 
that are independent of r = |x| outside some finite ball and for potentials with Cartesian 
anisotropy, i.e., for potentials for which \im X:j ^± 00 V(x.) exist for each j £ {l,...,n} 
Typically a scattering system containing a highly anisotropic potential involves a multi- 
channel structure. 

Time delay, for potentials vanishing at infinity, was first considered by Eisenbud jl^ . 
Bohm and Wigner By using asymptotic properties of the solutions of the stationary 
Schrodinger equation, they found that the energy derivative of the scattering phase shift 
may be interpreted as a time delay. Somewhat later Smith Q| suggested, as we mentioned 
at the beginning, to consider the excess sojourn time tx in a large spatial region X and 
to define time delay as the limit of Tx when this region tends to the entire configuration 
space M. n . He showed, also in a stationary framework, that this leads again to the Eisenbud- 
Wigner expression. Smith's proposal was formalized in the framework of time-dependent 
scattering theory by Jauch and Marchand Qj. These authors realized that the verification 
of the existence of the limit of tx as X — > W l represented a quite delicate mathematical 
problem. Later on a fair number of publications dealt with this problem; we refer to the 
review of Martin for details and references and mention that a satisfactory solution was 

2 [r^ . 
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In the present paper we shall consider the one- dimensional anisotropic case, i.e., Hamil- 
tonians of the form 

H = H + V (1) 

acting in the Hilbert space 7i = L 2 (R). Here Ho = P 2 = —d 2 /dx 2 is the usual free 
Hamiltonian and V = V{Q) is given by a real-valued potential V(x) assumed to have 
different limits at x — — oo and at x — +00; these limits will be denoted by Ve and V r 
respectively (see Fig. Here we have set h = 1 for Planck's constant and m = 1/2 for 
the mass of the particle, and we have written P and Q for the momentum and the position 
operator respectively in H. 
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FIG. 1: A typical potential V. 



Hamiltonians of this type present a two-channel structure and can serve as models in 
the theory of mesoscopic quantum systems. Scattering theory for such one-dimensional 
Hamiltonians, with potentials having different limits on the left and on the right, has been 
investigated mostly in the time-independent formalism [3], As regards the time- 

dependent approach, the existence of the M0ller wave operators was established in [21], and 
a complete and detailed mathematical study of the questions (i)-(iii) mentioned before can 
be found in a paper by Davies and Simon (2^]. As far as we know, time delay has been 
considered only for potentials for which Ve = V r , e.g., in j^j; of course this special situation 
is also covered by our results. 

We now describe the two representations of time delay in the above context, referring to 
Section for more details on our formalism. Because of the anisotropic structure of the 
potential there are two free Hamiltonians Hi = Hq + Ve and H r = Ho + V r entering into the 
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asymptotic description of the scattering states of H and hence also into the definition of S. 
If if) G Ti is a scattering state of H , given by a normalized vector and interpreted as the 
state of a particle at time t = 0, then its sojourn time (or dwell time) in the region [-R, R] 
(0 < R < oo) is given as follows: 

/oo pR 
dt dx \(e- iHt i;)(x)\ 2 . (2) 
-oo J-R 

Similarly one can introduce free sojourn times with respect to Hi and H r . More specifi- 
cally, let ip G 7i be the initial state corresponding to if>, i.e., satisfying ip = Q~<p, where Q~ is 
the M0ller wave operator for the limit t — > — oo. The incoming state tp may be decomposed 
into a part incident from the left and a part incident from the right: tp = tp%-\- tp r . By using 
a similar decomposition of the associated outgoing state Stp into a part (Stp)e propagating 
to the left and a part (Stp) r propagating to the right, one can then define the free incoming 
sojourn time T^(p) and the free outgoing sojourn time T^ ut (y9) associated to the initial state 
tp by 

/oo pR 
dt dx\[e- lH ^e + e- lHrt tp r ](x)\ 2 , (3) 
-co J-R 

/oo pR 
dt / dx |[e-^(SV), + e-^\S^)r]{x)\ 2 . (4) 
-co J-R 

The time delay in the interval [-R, R] induced by the presence of the scatterer (repre- 
sented by the potential V) is defined as the difference between the sojourn time of if> and 
the free sojourn times: 

7£*fo0 = T R (n-<p)-TF(<p). (6) 

As will be seen in Section these quantities are well defined for finite R. However for 
general states ip G Ti they are divergent as R — > oo, except when Vt — V r in which case they 
converge to the usual global time delay. The divergence when Vi ^ V r is not surprising: since 
the scattering is partially inelastic, the velocity in the state e~ lHt Q~ipe may be different for 
example from that in the free state e~ iHet ipi (assuming the particle is incoming from the left) 
at large positive times, so that the local time delays T^(tp) and r^ ut ((/?) will be proportional 
to R. A finite limit in this case, representing the global time delay for the initial state tp, is 
obtained by starting from the following symmetrized expression for the local time delay: 

tr(<p) = \ WM+^M] • (7) 
5 



It is interesting to know that, even in anisotropic cases (Vg ^ V r ), the (symmetrized) 
global time delay is identical with the time-independent Eisenbud-Wigner expression of 
time delay which is often employed in calculations and introduced by a somewhat formal 
argument using the notion of group velocity [jj. In particular, this identity permits one 
to have a time- dependent interpretation of the latter (see Section fVDjl . In terms of the 
S'-matrix S(E) at energy E, the Eisenbud-Wigner time delay operator at energy E is given 
as 

T(E) = -iS(E)*^ . (8) 

The family {T(E)} determines a self-adjoint operator T in the Hilbert space L 2 (M), and 
the Eisenbud-Wigner time delay is defined as the expectation of this observable T in the 
initial state tp: 

/oo 
(<p(E)\T(EME))dE, (9) 
-oo 

where (-|-) denotes the scalar product in L 2 (R), (-|-) that at energy E (details are given in 
Section ES| . 

To end this introduction we present our assumptions on V and outline the organization 
of this paper. We assume that the potential V is given by a real- valued function satisfying 

\V(x)-Vt\ < M(l + \x\)-" forx<0, (10) 
\V(x) - V r \ < M(l + for x > , (11) 

for some positive constants M and /i. Here Vi and V r are real numbers, with Vi < V r . The 
constant fi specifies the rate at which V(x) approaches its asymptotic limits. We assume 
throughout the paper that \i > 1 (short range condition). Some results will be derived only 
under further conditions on /i; the strongest hypothesis used is \i > 5 (Section IVBjl . 

The next three sections are devoted to time- dependent scattering theory. Rather than 



citing results from the paper by Davies and Simon [22|, who based their analysis on the Kato- 
Birman theory for trace class operators, we give a presentation in the framework of the more 
recently developed technique of differential inequalities (also called Mourre theory). This 
method, described in Section El will lead to various estimates on the rate of decay of wave 
packets at large times and then, in Section UTTI to the characterization of channel subspaces. 
In Section |IV] we discuss wave operators, asymptotic completeness and the S'-matrix S(E). 
Finally, in Section IS we use scattering theory to prove the existence of the limit defining 
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the global time delay and show its identity with the Eisenbud-Wigner representation. A few 
technical points will be explained in the Appendices. 



II. TIME DECAY OF WAVE PACKETS 

Time-dependent scattering theory is based on properties of the time evolution of wave 
packets and observables at large times t. One has to know that such quantities have limits as 
t — > ±00 or that they decay sufficiently rapidly in time (e.g., the integrals in ©-© defining 
sojourn times in bounded regions of configuration space should be finite). We describe here 
some basic estimates on time decay that will be used in the subsequent developments. The 
derivation of these results will be given in the Appendices. 

We write H K = H + k for the Hamiltonian with a constant potential given by the real 
number k and denote by T*p or tp the Fourier transform of a wave fonction p: 

{Fp){p) = 0(p) = ^=T e~^p{x) dx . (12) 

V ATX J-oo 

The spectral properties of H and H K are of course well known. The spectrum of H K is 
purely (absolutely) continuous and covers the interval [ac, +00). For the class of potentials 
considered here, i.e., satisfying (fTU | -([TT ]) with Vg < V r and /i > 1, the spectrum of H consists 
of an (absolutely) continuous part coinciding with the interval [Ve, +00) and possibly a 
set of non-degenerate eigenvalues < Vf, if /1 > 2, the number of eigenvalues is finite 



24j . We shall denote by 7i p {H) the subspace of the Hilbert space U = L 2 (R) spanned 
by the eigenvectors of the Hamiltonian H and by 7i c {H) its orthogonal complement in 
7i, corresponding to the continuous spectrum of H. The wave functions in 7i c (H) are 
the scattering states of the Hamiltonian H, also characterized by the property that they 
disappear at large positive and negative times from each bounded region in configuration 
space: 

ip E H C (H) <J=^ / \(e- lHt ^)(x)\ 2 dx ^ (13) 

as t — > ±00 for each fixed R G (0, 00). For a Hamiltonian H K with constant potential the 
subspace H C (H K ) is of course the entire space L 2 {R), and H P (H K ) consists only of the zero 
element of H, i.e., 7i p {H K ) = {0}. 
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A. Time decay in a constant potential 



In configuration space the time evolution of a wave packet ip in a constant potential is 
explicitly given as follows in terms of the free propagator (Lemma 3.12 in j^|): 



„—int poo /!/»._ „l 2 



iptix) = (e-^V)(^) = ^= J exp [ii—jf*-) <P(V) d y ■ ( 14 ) 

From this expression one easily obtains the following formula for the associated evolution in 
the Heisenberg picture of an observable f(Q), where / is a function of a real variable: 

e iHKt f(Q)e- iHKt = e- iQ2/4t f(2tP)e lQ2/4t . (15) 

As exp(ix 2 /At) converges to 1 when t — > ±oo, this equation expresses the fact that, at large 
times, position behaves approximately as momentum multiplied by t/m. 

bmce exp (iH K t) and exp(-iQ 2 /At) are unitary operators, one obtains from (|T5|) the fol- 
lowing important identity (by taking into account the Plancherel theorem, i.e., the unitarity 
of the Fourier transformation T): 

/oo ^oo 
\f(x)<p t (x)\ 2 dx= / \f(2tp)\ 2 \[^ Q2 ^ V )}(p)\ 2 dp. (16) 
-oo J — oo 

This can be used for obtaining time decay, in the L 2 -sense, of quantities of the form 
f(Q)e~ lHKt ip under assumptions on the momentum distribution of tp. 

We shall need the following estimates which can be obtained from (|16J) by taking for / 
a function of the Heaviside type or a function related to the potential V (see Appendix A). 
We say that tp is a wave packet with momentum in a subset A of the real line if <p(p) = for 
all p (jL A. We are particularly interested in decay properties at positive or negative times 
of wave functions with positive or negative momentum or with no non-zero momentum 
components close to p = 0. We have: 

(a) Let (p be a wave packet with positive momentum (i.e., with momentum in (0, oo)), Xo a 
real number and 6 > 0. Then there exists a constant Cg so that for t > 0: 

/xq poo 
\p t (x)\ 2 dx<C e (l + ty e / |(l + |x|)V(^ + ^o)| 2 ^ (17) 
-oo J — oo 



and for t < 0: 
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/OO 
|(l + |x|)V(^ + ^o)| 2 ^. (18) 
-oo 



(b) Let (p be a wave packet with negative momentum (i.e., (p{p) = for p > 0), Xq £ R and 
9 > 0. Then there exists a constant so that for t > 0: 



\<Pt{x)\ 2 dx < C e {l+ty 



•I'D 



and for t < 0: 



|^(x)| 2 c/x<C e (l + |t|)- 



|(1 + \x\) e cp(x + x )\ 2 dx 



|(1 + \x\) e ip(x + x )\ 2 dx . 



(19) 



(20) 



(c) Let / be a function satisfying |/(x)| < C(l + |a;|) _At for all ieR, some constant C and 
some /i > 0. Let ip be a wave packet with momentum in R \ (— Po,po) for some po > 0, and 
let 9 > 0. Then there exists a constant so that for all t £ R: 

1 1 



|/(x)^(x)|^x<C e 



|(l + |x|)V(z)rdx. (21) 



,(l+|t|) e (l + 2po|t|)2M_ 

Of course the above estimates are useful only if the integrals on the r.h.s. are finite. This 
requirement essentially amounts to a differentiability property of the Fourier transform of ip 
(<p should be 9 times differentiable in some sense). 



B. Time decay in a non-constant potential 

Results on time decay of wave packets in a non-constant potential are not so easy to 
obtain. The time evolution is now given by the unitary operators exp(—iHt), and it is 
useful to relate them to the Green's operator (H — z)' 1 , with z £ C \ R. For e > and 
E £ R we set 

6 {e) (H -E) = ^-[(H-E- te)- 1 - (H - E + te)- 1 } . (22) 

Then 

8 {e) {H -E) = f- / e lEt e- iHt -^ dt . (23) 

^ J -co 

If ip is a (square-integrable) wave function, we set ip t — e~ lHt ip and obtain from ()23|) by 
using the Plancherel theorem that, for any bounded operator B: 

/oo noo rcc poo 

dt dx e- 2eltl \(Bip t )(x)\ 2 = 2n dE dx \[B5 {e) (H - E)^](x)\ 2 . (24) 
•oo J —oo J —oo J —oo 

Let us consider wave functions ip which have non-zero components only in some finite closed 
interval A = [a, j3] in a representation where the Hamiltonian H is diagonal; such wave 



functions ip will be said to have energy support in A (with respect to H). For such a wave 
function and for energies E not in A, S^(H — E)ip converges to zero as e — ► (as S( e )(z) 
approximates the Dirac delta function), so that [2q 

/oo /*oo /*/? roo 

dt dx \(Bip t )(x)\ 2 = 2vrlim / ciE / c/x | [B5^ (H — E)ip] (x)\ 2 . (25) 
-oo J —oo e— >0 J a J —OO 

If one knows that there is a constant C such that 

/oo 
\[B5 (e) (H - E)^}(x)\ 2 dx < C (26) 
-oo 

for all energies E in A and all e in an interval (0, eo) for some eo > 0, then ()25)l implies that 
the L 2 -norm of Bip t is square-integrable over time, i.e., 

/oo /*oo poo 

dt dx \(Bi/j t )(x)\ 2 = / II^B^II 2 ^ < 2tt(/3 - a)C < oo . (27) 
-oo J —oo J —oo 

In (12 7J) and occasionally also further on we use the notation ||0|| for the Hilbert space norm 
of a wave function <f), i.e., \\<p\\ = [(^l^)] 1 ^ 2 = [j^ \<p{x)\ 2 dx] 1 / 2 . 

Bounds on L 2 -norms of the form (|26J) are a consequence of a Mourre estimate. Roughly 
the validity of a Mourre estimate means that, for some suitable self-adjoint operator A, the 
commutator between the Hamiltonian H and iA is strictly positive on wave functions ip of 
the type used above (i.e., having non-zero energy components only in A). More precisely, 
a Mourre estimate is satisfied on an interval A if there are a self-adjoint operator A, a 
compact operator K and a real number A > such that, for each ip having non-zero energy 
components only in A, the following inequality holds: 

M[H,iAW) > WW) + WW) . (28) 

In general A will be an unbounded operator, so that some additional technical conditions 
must be imposed in order for the commutator term occurring on the l.h.s. of (|2*Hj) to be well 
defined. 

Mourre theory, based on an inequality of the type (}2*Hj) . is an abstract method for studying 
general self-adjoint operators H. For the Hamiltonians considered here (i.e., one-dimensional 
Schrodinger operators) a suitable operator A is given by A = (PQ + QP)/4. As shown in 
Appendix B, a Mourre estimate is then satisfied on each interval above Vt, disjoint from the 
scattering thresholds, more precisely on each interval A = [a,(3\ with Vi < a < (3 < V r or 
V r < a < (3 < oo. 
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We mention some interesting general consequences of a Mourre estimate 
interval A = [ao, A)] m the interior of A (i.e., a < a® < A < (3), the Hamiltonian H 
has at most a finite number of eigenvalues (bound states), each at most finitely degenerate, 
(ii) The continuous spectrum of if in A is absolutely continuous (H has no Cantor-type 
spectrum in A), (iii) If Ao = [ao, A)] is an interval in the interior of A and disjoint from 
the eigenvalues of H, then the norm of the Green's operator (H — E — ie)" 1 , sandwiched 
between two operators (1 + |v4|) _1 , remains bounded near the real axis (i.e., for small values 
of e) at all energies E in A : 

sup ||(1 + \A\)~\H — E — iey\\ + lAir 1 !^) < C < oo , (29) 

EeA ,e^O 

where || • \\b(h) is the operator norm in the space B(Ti.) of bounded operators acting in 7i. 
A consequence of f!29|) is the validity of a strong version of the propagation estimate (|27|). 
with B = (1 + |^4|) _1 and wave packets ip having non-zero energy components only in Ao, 
namely (Proposition 7.1.1 in j^ ) 

/oo 
dx \(Bt/j t )(x)\ 2 < 8C \m 2 . (30) 
-oo 

In our situation the properties (i) and (ii) (even absence of bound states in (^,00)) 
were established independently of a Mourre estimate, as already mentioned before. However 
property (iii) is crucial for the approach to scattering theory presented below. Indeed, if (|2U|) 
is satisfied with the operator A indicated above, i.e., A = (PQ + QP)/4, then it also holds 
with (1 + l^l)" 1 replaced by (1 + |<5|) _1 (with possibly some different constant Co This 
implies the validity of the propagation estimate fpTOf) with B — (1 + |Q|) _1 . Consequently, 
if / is a function of a real variable satisfying |/(x)| < C(l + ) 1 for some constant C and 
all x, then there is a constant C\ such that for each with energy support in A : 



dt 



/oo 
\\f{Q)M 2 dt<Ci\\ii>\\ 2 . (31) 
-00 



In particular, for the Hamiltonians considered here, the time integral in (|3~T|) is finite for each 
scattering state xjj of H with bounded energy support disjoint from the scattering thresholds 
Vi and V r . 
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III. LARGE TIME LIMITS 



The bounds on certain time integrals obtained in the preceding section can be used to 
prove the existence of limits, as t — » ±00, of operators of the form e lhlt Je~ th2t , where J is 
a bounded operator and hi, h% are Hamiltonians. It will be enough to consider situations 
where J = g{Q) is multiplication in L 2 (R) by a function g(x). We are particularly interested 
in the cases g = 1, g — Xt an d g = Xn where Xi( x ) = 1 for 2 < 0, Xe( x ) = for 2 > 
and Xr — 1 — Xt- The functions xi an d Xt represent localization on the left and on the right 
respectively. For technical reasons we shall also introduce smooth approximations of xe an d 
Xr- We use the notation xt f° r the function xt an d for the operator xt(Q) °f multiplication 
by this function. 

All limits will be strong limits. We recall that, if {Wf}t e ]R and are bounded operators, 
then s-\im t _> +OQ W t = means that, for each wave packet tf>, W t ip — converges to 

zero in the Hilbert space norm, or equivalently that {W 7 ^} is Cauchy in the Hilbert space 
norm as t — >• +00. We first derive a useful formula (Eq. (}35|) ) for verifying that {Wt^} has 
this property. For < s < t we have in the Hilbert space norm: 



\W t ^-WM = sup 1(^1(^-^)1 
11011=1 

* d 



sup 

11*11=1 



— {(j>\WMdr 



(32) 

If W T = e ihlT Je- ih2T with h x = -d 2 /dx 2 + V 1 (x), h 2 = -d 2 /dx 2 + V 2 (x) and J = g(Q), then 



4-W T = ie ihlT (hiJ - Jh 2 )e- ih2T (33) 
dr 



and 



h x J - Jh 2 = -g" - 2g'^- + {V x - V 2 )g 

= -g"{Q) - 2ig'(Q)P + (Vi - V 2 )g(Q) . (34) 

Hence 

\\W t ^-WM< sup / \(e-^ T <f ) \i-g"(Q)-2zg'(Q)P+(V 1 -V 2 )g(Q)}e- lh2T ij}\dT . (35) 

||0||=1 Js 

The inequality ()35|) allows one to infer that {Wt^>} is Cauchy as t — > ±00 if a suitable 
estimate on the time decay of the integrand is available. If V 2 is a constant potential, (|33|) 
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corresponds to the well-known Cook method for proving the existence of limits of the type 
considered here. In this case it suffices to know the simple decay estimates of Section III Al 
an application concerns the existence of the M0ller wave operators in Section II V CI More 
refined techniques are needed for estimating the integral in ()35)) if V2 is a non-constant 
potential. Using estimates of the type (jSH), deduced from Mourre theory, we shall obtain 
two important results for scattering theory (existence of the limits in ()41|). (j^j) and (HH]))); 
details are presented in Appendix C. 

A. The channel subspaces 

In the simple case where h\ = h 2 = H K and J = xti convergence can be obtained without 
making use of an estimate of the form ()35|) . To know that e tHKt Xee~ tHKt converges strongly 
as t — > +00, it suffices to show that s-lim^ +00 e iiTK 'x£e _liTK V exists for a dense set of wave 
packets ip. We consider the following dense set: <p = ip + + <p_, where <p + has positive 
momentum, has negative momentum and (1 + \Q\) tp± are square- integrable for some 
9 > 0. The limit as t — > +00 of Xi e ~ * t( P+ is zero, which expresses the fact that e~ 
propagates towards the right (take Xq = and 9 > in (|17j)). Similarly ()19j) implies that 
s-lim i ^ +00 x r .e _i//K V- — 0. Since e %Htit is unitary, it follows that e iHKt xee~ tHKtl P- = L P- — 
e lHKt Xr e ~ V- converges to as t — > +00. In conclusion: s-lim^ +00 e* HKi x^e _;iHK V = </?_. 

As an operator, s-lim t _ >+00 e* HK 'x^e~* J?K * represents the observable of localization on the 
left at t = +00 in a constant potential; it does not depend on k and will be denoted by 
F^ t So is the (orthogonal) projection onto the subspace Ti.Q £ of wave functions that 
are localized on the left at t = +00 (in a constant potential), and it coincides with the 
projection II_ onto the subspace 7i_ of wave functions with negative momentum. 

One can similarly obtain the existence of the following limits and relate them to IT_ or 
to the projection Il + onto the subspace 7i + of wave packets with positive momentum: 

F + = s-lirn^ e iH ^x e e- iHKt = n_ , (36) 

F ~ = s-lim e iH ^x e e- iHKt = 11+ , (37) 

F+ r = ^lim e iH ^x r e- lHKt = Tl + , (38) 

F - = s-lim e iHKt Xr e- iHKt = n_ . (39) 

' t — ► —OO 
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In terms of these operators, the propagation properties of e lHRt tp± pointed out above may 
be expressed as follow: 

^lirn X ee- iHKt F+ r = ^lim ^e^F* = . (40) 
B. Subspaces of scattering states 

In the presence of a potential it is also possible to divide the set of scattering states 7i c (H) 
into two mutually orthogonal subspaces Ttf and 7i+ containing the state vectors localized 
on the left and on the right respectively at t = +00. For this one shows that the strong 
limits of e %m xi^~ %m and e %m Xr&~ %ut as t — > +00 exist on 7i c (H) and define two projections 
F^~ and F r + in TC C (H), with F^F+ = 0. A similar decomposition TC C (H) = TiJ © TC~ 
exists, corresponding to the limit t — > — 00. In contrast to the case of a constant potential 
(where T~Cq £ = / Hq 7 ), these two decompositions of 7i c (H) are different in general; indeed the 
equality of TiJ and 7i+ would mean that the potential is reflectionless. On the other hand 
and Hj (and similarly 7i+ and 7i~) are related by time reversal. The antiunitary time 
reversal operator 0, given by complex conjugation, i.e., (®ip)(x) = ip(x), commutes with 
the Hamiltonian H and with xii so that Qe lHt XeG~ tHt ip = e~ lHt XeG lHt Qifj. It follows that i/j 
belongs to 7if if and only if ©■?/> belongs to H.J . Obviously these relations also hold for Ti^g 
and T~C^ r , i.e., QTi.^ e = Hq £ and QH^ r = H^ r . (An equation bearing double signs always 
has to be interpreted as two independent equations, one for the upper and one for the lower 
sign). 

Denoting the projection onto the subspace TC C (H) of scattering states by F C (H), we set 

Ff = s-lim e iHt X£ e- iHt F c (H) , (41) 

t — ► ± OO 

F± = s-lim e iHt x r e- iHt F c (H) . (42) 

t — > i OO 

The existence of these limits can be obtained by using dHHJ) with W T = e iHT g(Q)e- iHT , where 
g is a smooth approximation of xe. or Xr- Details on this are given in Appendix C. Below 
we deduce alternative expressions for these limits (Eqs. (|43j) and (|4l)j) ) and determine their 
properties. We shall consider F^~. 

Let F P (H) = J2k=i lv 9 fc)(v 9 fc| be the projection onto the subspace 7i p {H) of bound states 
of H\ (fx, . . . , (pw are normalized eigenfunctions if H has N bound states (we consider the 



14 



case where N ^ 0). Denoting the eigenvalues by E±, ... , we have e * H Vfc = e lEkt{ Pk-, 
and the square of the norm of F p (H)e tHt Xe^~ lHt ' l f ; is just 

TV N 

£|<^|e^e-^V>| 2 = EK^IX^W • ( 43 ) 
fc=i fc=i 

We set ip t = e~ lHt ip and define the function X{a,p] by X(a,/3]( x ) = lifa<x<[3 and 
X(a,/3](^) = otherwise. Then \i = X(-oo,o] = X(-oo,-R] + X(-R,o] for any R > 0, hence 

(Pklxe^t) = (X(-oo + MX(-R,o]^i) • (44) 

By the Schwarz inequality, the absolute square of the first term on the r.h.s. of 
is bounded by J_^\(pk(x)\ 2 dx \ifj(y)\ 2 dy, which can be made arbitrarily small by 
choosing R large enough. The absolute square of the second term is bounded by 
JToo Wki.x)^ dx f® R \i/j t (y)\ 2 dy, which converges to zero as t — * +oo by (JT3J) if ip G H C (H). If 
N is finite, we conclude that F p (H)e iHt x e e- iHt F c (H) converges strongly to zero as t — > +oo 
[29^ . Since F P (H) + F C (H) = 1, it follows that the limit defining Ff exists if and only if 
s-\im t ^ +oc F c (H)e iHt x e e~ iHt F c (H) exists, and that 

F+ = F C (H)F+ = t s-liin F c (H)e iHt X ee~ lHt F c (H) . (45) 

The step function xi m Q45j) may be replaced by a smooth approximation g. Let g be a 
smooth function satisfying < g{x) < 1 for all x, g(x) = 1 for x < — 1 and g{x) = for 
x > 0. Then, by (|13|) with i? = 1, x& e h> ~ g{Q)e~ lHt ip converges to zero as t — > +00 for 
each ip e H C (H), hence 

F+ = s-lim F c (H)e iHt g(Q)e~ iHt F C (H) . (46) 

£ — > +00 

Let us mention some simple consequences of (|43j) : (i) is self-adjoint, i.e., (F^)* = 
(ii) F £ + is idempotent. Indeed, using (i), the unitarity of e lHt , the fact that xj = Xi an d the 
identity (jlSJ), one obtains 



M(i?)V> = {F?ip\F?1>) 



—iHtp ^ff^jtp^ 
JHt„ , -iHt 



\im{e im xee- im F c (H)ip\e iHt Xe e- iHt F c (H)^) 

t— »oo 



lim( V p|F c (^)e^ t x^^^ c (^)^> 



t— »oo 



<^|F+^) • (47) 
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Hence {F^) 2 = . The properties (i) and (ii) mean that F^~ is an orthogonal projection in 
H C (H). The subspace onto which it projects is denoted by TC^. This subspace is invariant 
under the evolution: e~ tHs leaves invariant (s G 1), or equivalently F^~ commutes with 
e~ lHs (and hence with H). Indeed, since e~ lHs commutes with F C (H): 

F+e~ iHs = s-lim e iHt x e e- iH{t+s) F c (H) 

t — ► + OO 

= s-lim e iH{T - s) X£ e~ iHT F c (H) 

T — » +OO 

= e~ iHs F+ . (48) 

Also, by an argument as in (J4T)) . using the relation \rXt = 0; one finds that F+F^~ = or 

s-lim Xr e- iUt Ff = s-lim Xe e~ iHt F r + = . (49) 
t — > +00 t — » +00 

This shows that is orthogonal to H/, and clearly 7i c (H) = © Tij since F^~ + F r + = 

IV. SCATTERING THEORY 

As explained in the Introduction one needs to answer some dynamical questions (points 
(i) to (hi)) in the framework of scattering theory in order to define rigorously the notion 
of time delay. We first introduce the asymptotic condition representing the fundamental 
idea of scattering theory. This leads us naturally to the introduction of the M0ller wave 
operators Q and then the scattering operator S. The anisotropic structure of the systems 
considered here is reflected in the appearance of channel operators. Existence, properties 
and the interpretation of these channel operators are then discussed. Finally we introduce 
the S"-matrix S(E) entering the Eisenbud-Wigner representation of time delay. 

A. Asymptotic condition 

The fundamental idea of scattering theory, expressed in the context of anisotropic systems 
considered in this paper, is that at large (negative and positive) times t a particle in a 
scattering state ip t = e~ tHt tp G 7i c (H) is located in a region far from the scatterers, where 
the potential V is essentially constant (approaching Vi at x = —00 and V r at x — +00), 
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and therefore should behave essentially as a free particle (evolving with Hi = H + Vt and 
H r — Hq + V r respectively). 

This idea is called the asymptotic condition and is formalized as follows. Let ip E TC C (H) 
be a scattering state of H. Then there should exist free scattering states iff E Ti. of Hi and 
tpf ETi of H r so that 

lim / \(e- iHt ^)(x) - (e- lHet <pf)(x) - {e- iHrt <pf)(x)\ 2 dx = , (50) 

where as before this condition has to be interpreted as two independent relations, one for 
each sign in ±. 

B. Wave and scattering operators 

To show that the conditions (JHOJ) are satisfied, we introduce the Mfiller wave operators 

o± = nf + nf (51) 

with 

flf = e iHt e- iHet F± e , (52) 

ft± = s-lim e iHt e- iHrt F± r , (53) 

t — ► i oo 3 

where the projections F^ e and F^ r are those defined in (j3T)j) - (f3T?j) . We call the operators Clf 
and Qf the channel wave operators. 

It is straightforward to see that the asymptotic condition (J50|) is equivalent to the exis- 
tence of the wave operators Q^ 1 and their completeness, i.e., their ranges T^f^) are the entire 
continuous subspace TC C (H), justifying a posteriori the introduction of the projections F^ t 
and F^ r selecting the states in Ti which have the appropriate asymptotic behavior. These 
facts will be proved in the next subsection. Assuming this one can define the scattering 
operator by 

S = . (54) 

The decomposition (JSTj) of the wave operators into a left and a right part leads to the 
following structure of S: 

S = Si + S r = Su + Sri + S rr + Sir , (55) 
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with Si = Su + S r £, S r = S rr + S^ and S a b = where a and b stand for £ or r. We 

call the operators Sat, the channel scattering operators. Their properties and interpretation 
are established in the next subsection. In the sequel we shall use the subscripts a, b, c and 
d to denote £ or r. 

The introduction of the sum of the channel wave operators, Eq. (j51j) . may seem somewhat 
unusual for a multichannel system. In our context this sum represents an interesting and 
useful operator, as shown under point (iii) of the next subsection. It should be noted 
that the channel structure arising here is different from that occuring for example in the 
quantum-mechanical iV-body problem. In the latter the channel subspaces may overlap but 
are independent of the sign of time, whereas for the scattering systems considered here the 
channel subspaces depend on the sign of time (see (jHEJ-fjUHl)) but represent, for each sign 
of time, a decomposition of the Hilbert space 7i into mutually orthogonal subspaces. The 
wave operators Or in (|51|) may be expressed in terms of the free Hamiltonians H m and H ont 
introduced in (}65|) as follows: tt~ = s-\im t ^-ooe im e~' tHint and f2 + = s-\im t ^ +00 e lHt e~ %H °^ t . 
[The quantum-mechanical iV-body problem is often written in a similar form, called the 

n n 

two-Hilbert space formulation (see e.g., [21] or by introducing an auxiliary asymptotic 
Hilbert space 7Y as (the orthogonal direct sum of the channel subspaces) and an asymptotic 
free Hamiltonian H as acting in 7i as and independent of the sign of time.] 



C. Existence and properties of the channel wave and scattering operators 

In this subsection we show that, if one assumes the potential V to satisfy (fHH) - (fTTf with 
/Lt > 1, then the channel wave operators are well defined, i.e., the time limits involved in 
their definition (f52|) - ([53|) exist. We then establish their properties and those of the channel 
scattering operators S a b (assuming that \i > 2). 

The existence of the channel wave operators can be obtained by invoking only properties 
of time decay in a constant potential. We show the existence of f^j", the other channel wave 
operators can be handled similarly. Let Q(t) = e lHt e~ lHet F^~ e and write 

Q(t) = e lHt X ee- iHit F+ e + e iHt Xr e- iHet F+ e , (56) 

where xe an d Xr represent the configuration space localization on the left and on the right 
respectively introduced in Section UTTl 
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Using the unitarity of e lHt and (|40|) it follows that the second term on the r.h.s. of (|56|) 
vanishes as t — > +00, so we have 

fi+ = s-lim e im x^e- iH£ 'F + 



t 

= ^s-lirn^ e iHt g{Q)e- mt F+ e , (57) 

where g is the smooth approximation of \i introduced before Eq. (|46j). 

We consider the following dense set of wave packets ip in the subspace H^f <p has 
(negative) momentum in a bounded closed set not containing p = and is such that (p 
is three times continuously differentiable. From (|H5jl with h\ = H and /12 = -f^, hence 
W T = e iHr g(Q)e- iHtT , one obtains that 

||Wip- < / iV T rfr , (58) 

where 

iv T = ||[-/(g)-2V(g)p + (T/-^(g)] e -^>|| . (59) 

Since g' and g" vanish outside the interval [—1, 0] and |[V(a;) — Vg] ^(o:) | < M(l + |x|)~ M with 
/i > 1, the integral J °° iV r <ir is finite as a consequence of (J2T|) (with 6* = 3). Thus exists. 

Next let W(t) = e lHet e~ lHt F^~ . As shown in Appendix C, these operators are strongly 
convergent as t — > +oo. By arguing as in (|5fjjl - (|57jl and using ()49jl one obtains the following 
expression for their limit: 

W = s-lim w(t) = s-lim e iHit g(Q)e- iHt F+ . (60) 

t — ♦ + oo t-^ + oo 

Let us show that the operator W is in fact the adjoint (fi/)* of the channel wave operator 
ilf. Using the expressions (JH7j) and (160( 1 for f2^~ and W respectively one obtains the following 
useful identities by proceeding as in ([47)1: F^~Q~^ = and = W. Then 

((sv)vi^> = (mtv) 
= (F e + mi<p) 

= lim (F Q y H ^e- im F^\ V ) 

= (wvb) • 

Hence (ft/)* = W. 
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We now collect the properties of the wave and scattering operators. 

(i) Cl^f is isometric on Ti.Q a . Indeed, we have H^V^II = ^ m t^±oo \\ elHte ~ lHat F^ a tp\\, which 
is equal to ||<^|| if (p G 7i^ a and vanishes if ip is orthogonal to H^ a . 

(ii) The range TZ(Q^) of fl^ is the whole subspace TL^ ■ Indeed, let us consider for example 
the case of f^j~. Notice first that the identity F^Q,^ = implies that 7^.(0^") is a subspace 
of Til . Then let if> e TC^ and set <p = Wip. From the identity F^ e W = W one deduces that 
ip belongs to H.q £ . Finally we have ip = VL^pi since 

= lim ||e^e-^F >-^|| 
= lim \\<p - e iHlt e- iHt if}\\ 

t— *+oo 

= \\ip-Wip\\ = o . 

(iii) The wave operators are isometries from 7i onto 7i c (H) and therefore are complete, 
i.e., 1Z(H ± ) = 7i c {H), where TZ(Q, ± ) denotes the range of Q ± . Indeed, writing f2 ± = £lf + Qf, 
using the points (i) and (ii), and recalling the decompositions H = 7Y^©7i^ r and H C (H) = 
7~if © 7~Cf established in Section IHII it follows immediately that fi 1 * 1 is an isometry mapping 
H onto H C (H). 

(iv) The channel wave operators satisfy intertwining relations: e~ lHt Q^f = Q^fe~ lHat . More 
generally such relations hold for a large class of fonctions <3>: = Q^Q(H a ). These 
relations express the conservation of energy (kinetic plus potential) in scattering processes 
(see the Remark in the next subsection). 

(v) S is unitary. Indeed, by (iii) we have 

s*s = (n-yn + (n + yn~ 

= {q~)*f c {h)q~ = (ir)*ir = / 

and similarly SS* = I, where / denotes the identity operator in Ti. 

(vi) From (i) one easily sees that the operators S a and S a b are not isometries. However, 
it is clear that S a is isometric as an operator from a to Ti^ © Hq t = 7i and that the 
channel scattering operators S a b associate to each initial state in b a final state in Ti.Q a . 
For example the operator Si maps initial states with positive momentum (corresponding to 
particles that are incident from the left) to the associated final states. If tp is a wave packet 
with positive momentum, then Slip is decomposed into the superposition of S^cp (the part 
of the associated final state propagating towards the left) and S r i<p (the part propagating 
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towards the right). The meaning of the operator S r is similar when applied to particles 
incident from the right. 

(vii) The channel scattering operators satisfy the following important identities: 



(s ab ys cd = o, s ba (s dc y = o (ei) 

if a 7^ c, and 

(S a b)* S ac = SbcF b, Sba(S ca )* = &bcF^ b ■ (62) 

a=£,r a=£,r 

Indeed, (JHH) follow directly from the fact that fi+(fi+)* = 5 ac F+ and fi~(fi~)* = S ac F~ . To 
check for example the first equation in (|62[) we write 

(S a b)*S ac = (Sib)* Sic + (Srb)* S rc 

a=£,r 

= (^ynKntyn; + (^yn+(n+yn; 
= (^)*f c (h)q: = (n-yn: = 6 bc F- b . 

(viii) The intertwining relations for the channel wave operators immediately transcribe 
into the following intertwining relations for the channel scattering operators: e~ tHat S a b = 
S ab e- lHbt . 

(ix) Let be the time reversal operator given by complex conjugation: (Oip)(x) = 
J{x). It satisfies 0F O ± 6 = F* a , 0F±0 = (see Section EH) and ee iHt e~ iHat F^ a 6 = 
e -iHt e iH a t F ^ hence e^±e = The last relation also implies that 6(^)*9 = (fif)*. It 
follows that 

es ab e = (s 6a )* , (63) 

hence, recalling the decomposition S = J2 a b=e r Sab, one obtains 

050 = S* . (64) 

D. Scattering matrix 

As observed in Section ITTT1 (see in particular (|3T|) and (jSHJ) the subspace 7i + of wave 
functions with positive momentum must be associated with different physical situations 
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when one considers incoming or outgoing states. The same is true for the subspace of 
wave functions with negative momentum. Wave functions in H + (or 7i_) describe states that 
are incoming from the left (or the right), their evolution being governed by the Hamiltonian 
Hi (or H r ). The component in 7i + (or Ti-) of an outgoing state corresponds to a particle 
propagating towards the right (or the left), its evolution at large times being governed by 
H r (or Hi). Thus it is natural to introduce two free Hamiltonian H m and H ont by 

H' m = HtIL+ + H r U_, H out = H r Tl + + #T1_ (65) 

(n + and Il_ being the orthogonal projections in Tt onto 7i + and H- respectively). 

To diagonalize H m we identify TC = H + © with a subspace TC m of the complex 
Hilbert space L 2 ((Vi, oo); C 2 ) of square-integrable 2-component functions of the variable 
E G (Vg, oo). The elements of H m have only one non-zero component for E G (Ve,V r ) 
and in general two non-zero components for E > V r , their definition being as follows. Let 
= 0+ + 0_ be a wave function (0 + G H + , 0_ G TiJ). We denote by (£> m (E) G C 2 the value 
of at E in 7i m . In consideration of the meaning of <j)± when is viewed as an incoming 
wave packet, we use the notations ^{E) for the component of <j) m (E) associated to 0+ and 
(j)f'(E) for that associated to 0_. So 



(E) = "< ~ (66) 




with 



t £ ) = p^jpS^) (^>^) (67) 
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C(£) = _ 1/4 0- (- ^ ~ K) (£>V r ) (68) 

and 0J. n (-EO = for < K- The normalization factors on the r.h.s. of (|67jl and (Jfi%|) are 
chosen such that the identification of H with the subspace 7i in of L 2 ((Vi, oo); C 2 ) is unitary. 
Clearly (H in (f)j m (E) = E(f m (E). 

The diagonalization of H ont is achieved similarly by identifying Ti with a subspace T~C out 
of L 2 ((V^, oo); C 2 ) in the following manner: for = + + 0_ interpreted as an outgoing state 
(0+ G H + , 0_ G H-) we set 

^T\E) = { ^ E l Vi)Y /M -VE^Vi) (E>V e ) (70) 
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and 

with (p° r nt (E) =Qj£E<V r . 

From the relations S a bH b = H a S a b one finds that 

SH in = H out S . (72) 

If S is viewed as an operator from 7i m to 7i out , (|72jl means that S maps the value at energy 
E of an incoming state to the value of the associated outgoing state at the same energy E: 
for each E > V r there is a (unitary) 2x2 matrix S(E), called the scattering matrix, such 
that 

(Sp) out (E) = S(E)p in (E) (73) 

or more explicitly 



(S<p)?*(m = (Srt(E) S rr (E)\ Uf{E) 
(Stp)r\E) " \S U (E) S er (E) U n (£) 



(74) 



For Ve < E < V r , where ^(E) = (Sip)° ut (E) = 0, there is a complex number Su(E) 
of modulus 1 such that (Sip)^ lt (E) = Su(E)<f 1 l l (E). The relation (|T3|) specifies the wave 
function Sip, viewed as an element of 7i out , in terms of the wave function ip represented in 
U m . 

We point out a simple way of arriving at the above structure. Consider the operator J 
in H that interchanges TC + and H- (i.e., the parity operator): 

W)(p)=k-P)- (75) 

It satisfies J* = J, J 2 = I, JH± = H T , and it interwines H in and H ont : 

H' m J = JH out . (76) 

It follows that H' m JS = JH ont S = JSH' m , i.e., JS commutes with H' m . So (see e.g., 
Proposition 5.27 in j^]), in the representation 7i m diagonalizing the self-adjoint operator 
H m , the unitary operator JS is decomposable, i.e., for example for each E > V r there is a 
unitary 2x2 matrix cr(E) such that for each (p in Ti: 

{JSpf\E) = a{E)^{E) . (77) 
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It is easy to check that, if written as an operator j7in->out from H m to H out , J is given by 

= (J^^rHE) = ( ° 1 I 4>'"(E) . (78) 

Hence, writing the scattering operator S : H m —> H out as S = J(JS) = Jm->o\it(<JS)i n ^- m , 
one sees that (for E > V r ) S has the form (|74|). with S'(.E') = ^(-E 1 ) a unitary 2x2 

V V 

matrix. 

We add a few comments concerning the S"-matrix S(E) in ()74j1 . For E > V r the complex 
numbers Sa(E) and S rr (E) are the reflection amplitudes at energy E, whereas S r i(E) and 
Se r (E) represent the transmission amplitudes. Time reversal symmetry and unitarity lead to 
important relations between these quantities. Equation ()64j) may be rewritten as S = 05*0. 
As an operator from 7i out to 7i m , S* is decomposable, its component at energy E > V r is 
just 



S r AE) Su(E)\ 

S*{E) = S{E)* = V ; V ; . (79) 
\S rr (E) S ir {E) I 

The time reversal operator 0, expressed as an operator from 7i m to 7^ out , is also decompos- 
able (C denotes complex comjugation): 



C \ I d) in (E) 

(00) out (£) = I 4> in (E) 



C OJ \<Pf( E ) y 
So, writing the operator QS*Q : TC m —> H out as 0i n ->out(£*)oiiWii0in->out> the equation 
S = 05*0 implies that 

S r e(E) S rr (E) \ I Se r (E) S rr (E) 
S ee (E) S er (E)J \S ee (E) S re (E) 

hence 

S r e(E) = S ir (E) . (80) 

The transmission amplitudes from left to right and from right to left are equal. The unitarity 
of S(E) then leads to 

\S ie (E)\ 2 = l-\S re (E)\ 2 = \S rr (E)\ 2 , (81) 
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i.e., the reflection probabilities at energy E > V r from the left and from the right are the 
same. 

Remark. The variable E on the r.h.s. of (fTH|) refers to the value of H m whereas that 
on the l.h.s. relates to H ont . As an illustration consider a scattering state ip in TiJ (i.e., 
incoming from the left) having energy support (with respect to H) in a very small interval 
A centered at Eq > V r . We shall say that tp energy Eq. By the intertwining 

relation (Qj)*Q(H) = &(Hg)(Qj)* the associated initial state if = in T~Cq£ is a state 

at free energy E (relative to H m , i.e., (p m (E) = fg n (E) = for E outside A); for this state 
it is natural to view Eq as being composed of a potential energy Vg and a kinetic energy 
Xg pd Eq — Vg. Likewise the components Sggip = fl^ip and S r £tp = of the final state 
S t (p are at free energy E (relative to H ont ), i.e., {Sggip) out (E) = and (S r£ ip) out (E) = for 
E ^ A, their kinetic energy being Xg w E — Vg and A r ps E — V r respectively. This expresses 
the fact that reflection is an elastic process while transmission is inelastic if Vg ^ V r (the 
operator Sgg commutes with H = P 2 whereas S r g does not if Vg ^ V r ). 

V. TIME DELAY 

As explained in the Introduction time delay expresses the excess time that scattered 
particles spend in the scattering region when compared to free particles and therefore is 
naturally formalized, in the framework of time-dependent scattering theory, in terms of 
sojourn times. However one often uses in calculations the Eisenbud-Wigner expression 
of time delay given in terms of the S'-matrix. In this section we prove that these two 
representations of time delay are in fact identical. We first show that for appropriate initial 
states ip the sojourn times (jSJ)-© are finite for each finite R, so that the local time delays 
(JSJ-Q are well defined. We then show that in the case Vg ^ V r only the symmetrized 
expression Q of local time delay has a finite limit as R — * oo; we call this limit the 
global time delay. In the context of time-independent scattering theory we first observe that 
the Eisenbud-Wigner operator T = {T(E)} is well defined and self-adjoint, so that the 
Eisenbud-Wigner expression of time delay Q is quantum-mechanically natural. We finally 
establish the main result of this paper, namely the identity between the global time delay 
and the Eisenbud-Wigner expression of time delay. 

We assume that fi > 2 in the decay assumptions (fTTijl - fTTf on the potential. In order 
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to avoid longer expressions, we discuss the case of wave packets that are incident from the 
left (<p G = H.+). The same type of arguments are applicable to wave packets incident 
from the right (ip G 7io r = H-) and to general initial states (ip G Ti). More specifically we 
shall consider initial states ip G 7i+ having energy support (with respect to He) away from 
the thresholds Ve and V r , and some decay in configuration space. For this we introduce a 
parameter 9 > (its values will be specified further on) and denote by T> 1 ^ the set of wave 
functions <p G 7i + satisfying 

(i) there are intervals Ai = [E\, E[] and A2 = [E2, E'^\ (depending on ip) with Vi < E\ < 
E[ < V r < Ei < E' 2 < 00 such that ip has energy support (with respect to He) in 
AiUA 2 , 

(ii) the following integrability condition: 



To handle the final states S^ip (which evolve with H ont ) we introduce a similar set of wave 
functions Vg ut as follows: = (f> + + 0_ G H + © 7Y„ belongs to Vg nt if there are intervals 
Ai and A2 as in (i) above such that: (1) <p + has energy support with respect to H r in 
A2 and (fi- has energy support with respect to He contained in Ai U A2, (2) + and 0_ 
satisfy the integrability condition (|82[). Notice that, if tp G 7i+ has the support property (i) 
in the definition of D™, then = Sgip satisfies the support condition in T>Q Ut (this follows 
from the intertwining relations H a S a b = S a bH b , observing that (f> + = S r i<p and 0_ = See^)- 
Nevertheless some assumptions on the decay of V are needed so that Si maps into Vg ut 
(see Section IV B|) . 

A. Sojourn times 

For < R < 00 the sojourn times in the interval [-R, R] associated to a wave function 
ip eH + are defined by (see ©-©) 




(82) 




-R 
r R 



(83) 



(84) 



j —00 



/ dt 



-R 



dx \(e 



r ^)(x) + (e-^5 r ^)(x)| 2 , 



(85) 
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If in addition tp has energy support (with respect to Hg) away from Vg and V r these quantities 
are finite. This is not surprising since such wave functions describe states of a particle with 
non-zero velocity. Mathematically the finiteness of T R (Q~ip) follows from (pTTf with f(x) = 1 
or if \x | < R or \x\ > R respectively (ip = VL~ip has energy support with respect to H away 
from Vg and V r by the intertwining relation HQj = VLJ Hg, and < (-R+l)(l + |x|) _1 ). 

The finiteness of T R n ((p) and T^ ut (y?) can be obtained in the same way, applying ()31|) for the 
Hamiltonians Hg and H r , with ip t = e~ tHet (p, ipt — e~ %Ult Sggty and ip t = e~ tHrt S r g(p. 

If f G TC+ is a wave packet satisfying (|%2~)) for some 9 > 1, the finiteness of T R (<f) is also 
an immediate consequence of the decay estimates (J 17)1 and (j!8j) (take xq = R in (fTTjl and 
Xo = — -R in (jl8|) ). Similarly the finiteness of T^ ut (yj) then follows from fTTj) - (f2T)j) provided 
that one knows that Sggip and S r g<p also have the decay property ()82)1 with 9 > 1. 

B. Local time delay 

The results of the preceding subsection imply that, for initial states belonging to T> 1 q with 
9 > 0, the local time delays ©-© are finite for each finite R. We now turn to the question 
of existence of a limit of these quantities as R — > oo. Following ^] we proceed in two steps: 
(1) approximate T^((p) and T^ ut (</?) by expressions giving the same limit but involving the 
scattering operator Sg rather than the wave operator Qj, (2) use an asymptotic expansion 
(for large R) of J ±o ° e tHKt X{-R,R){Q)e~ lHKt dt, where X(-r,r)( x ) = 1 or if |x| < i? or \x\ > R 
respectively. We treat here step (1) and discuss step (2) in the next subsection. We assume 
that (f belongs to T>™ for some 9 > 4. 

Let us consider r R n ((p). Setting tp t = e~ lHet ip we have 

/oo pR poo 

dt / dx [|(^)(*)| 2 - \ft(x)\ 2 ] = / I R (t) dt , (86) 
-oo J—R J —oo 

where /ij(t) denotes the integral over the variable x. Observe that for any finite to and ti 
(t < ti), the quantities dt j^ R dx\{Vt^ (f t ){x)\ 2 and dt J^ R dx |<^(x)| 2 are increasing 
functions of R each of which converges to (ti — 1 ) JZ° Iv 9 ^)! 2 dx as R — > oo. Hence, for any 
to < h- 

/to foo 
I R (t) dt+ I R (t) dt + o(l) as R -> oo . (87) 
-oo J t± 

Now, by the definition of the wave operators, f^V 9 * an< ^ ft approach each other (in the 
Hilbert space norm) as t — > — oo. One thus expects that the integral over (— oo, to) m fjHTj) is 
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negligible if to is negative and sufficiently large, so that it suffices to consider the contribution 
to Tr(<p) coming from large positive times. As explained below this is indeed the case. 

Denoting by X(-R,R) multiplication by X(-r,r)(x), writing in terms of inner products 
and then using the Cauchy-Schwarz inequality we have 

\I R (t)\ = \(SlJip t \x(-R,R)(£lJipt- <Pt)) + (ttjipt- <Pt\X(-R,R) l Pt)\ 

< W^IftW ■ \\X{-R,R)(^e<ft - <Pt)\\ + K<Pt - <Pt\\ ■ \\X(-R,R)<Pt\\ 

< 2\\(p\\ ■ \\(fij - I)<pt\\ . (88) 

We now observe that 

r° d 

(fi7 - I)ip t = - — \e iHu e- iHtU ] e~ iHtt ipdu 
J-oo du 

= -ie~ lHt [ e iHs {V -Vt)e- iHiS ^ds . (89) 



By taking into account the (continuous) triangle inequality in the Hilbert space norm, one 
obtains that 

< f \\(V -V t )<p a \\ds 

J — oo 

t pt 



< \\{V-Ve)xMs\\ds+ \\{V -V t )xrV.\\ds . (90) 

J — oo J — oo 

The norms appearing in the integrands can be estimated (for s < 0) by using (|21[) and (|18j): 

\\(v-v t )xe<p.\\ < c(i + \ s \yp 

\\(V-Vt) X r<P.\\ < C{l + \s\r e l 2 . 

Here C is some constant (depending on ip and 9) and p = min{/i, 6/2}. It then follows from 
(JHEj) and (jnOJ) that, for some constant C and t < 0: 

\I R (t)\ ^Clt] 1 -?. (91) 

This bound is valid for all R > (the constant C depends on tp and 9 but is independent 
of R). Since we assumed that 9 > 4 and /i > 2, we have p > 2, so that | J*^ J R (t)<it| can be 
made arbitrarily small, independently of R, by choosing t sufficiently large (negative). 

As t — > +oo, 1^7 V 9 * approaches S^t. Writing Qjip t — Se<pt = {^f— I)Su<ft+(^t ~ I)Sre<fit 
and proceeding as above (using the decay estimates fTTJ) . (JTJJ) and (J2H)) one finds that 



dt / dx[\(nj<p t )(x)\ 2 -\(s m )(x)\ 2 ] 

-R 
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(92) 



can be made arbitrarily small, independently of R, by choosing ti > large enough, under 
the proviso that Sgtp belongs to Vg ut for some 9 > 4. Hence, setting 

poo pR 

o%{<p) = dt dx[\(S m )(x)\ 2 -\<p t (x)\ 2 ] (93) 

JO J -R 

we have 

lim i%(<p) = lim o%(<p) . (94) 

R— >oo R-^oo 

The relation (|94|) shows that, as i? — > oo, either T^(ip) and ^(y?) converge to the same 
finite limit or both are diverging. 

By the same type of arguments one finds that 

lim r° ut M = lim oT{<p) , (95) 

R — >oo R— «X) 

with 

<(<f) = - f dtf ' dx [\{S m ){x)\ 2 - |^(x)| 2 ] . (96) 

The relations (|9*4^) and ()95|) are satisfied for wave functions tp belonging to X>™ for some 
> 4 such that S^(/3 belongs to Vg nt (for some possibly different value of 9 > 4). This 
restriction essentially amounts to a differentiability condition on the S'-matrix S(E), namely 
that S(E) should be 9 times differentiable (with respect to E) away from the thresholds 
Vt and V r . Differentiability of S(E) can be obtained for example from differentiability 
assumptions on the potential V or from assumptions on the decay of V at large \x\ (i.e., 
assumptions on the parameter \x occuring in (|10 p -([ll| ) ). The entries of the matrix S{E) 
are simple expressions in terms of the Jost solutions of the stationary Schrodinger equation, 
and it suffices to know that the Jost solutions are 9 times differentiable with respect to 
the energy parameter. We refer to Appendix D for a brief discussion of these questions. 
The simplest situation is obtained by taking 9 = 5. Then the above-mentioned conditions 
on (f and S^ip are satisfied if ip fulfills the support condition (i) in the definition of X^ n 
and (p is five times continuously differentiable, and if fi > 6 in ([10j ) -([li p (then S(E) is five 
times continuously differentiable away from the thresholds). If one admits differentiability 
of fractional order (Lipschitz or Holder continuity of ip and Sip), then it suffices to assume 
that [i > 5. If /j, > 5, then the fourth derivative S^(E) of S(E) is locally Holder continuous 
with exponent 7 = min{l,/x - 5}, i.e., \S${E) - (E')\ < C A \E - E'\i for E, E' in any 
closed energy interval A not containing the thresholds. If ip^ is also Holder continuous with 
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the same exponent 7, then for 5 < 7: 

/oo 
1(1 + \x\) s (l + \x\)\S ab p)(x)\ 2 dx < 00 , (97) 
-00 

i.e., Sftp £ The finiteness of the integral in (|97|) follows from a result in classical 

Fourier analysis stating that, if / £ L 2 (M), f(p) = for all p outside some bounded subset 
of K and \f(p) - f(p')\ < C\p -p'\~< for some < 7 < 1, then |(1 + \x\) s f{x)\ 2 dx < 00 
for 6 < 7 (see Section 4.13 in |31|, in particular Theorem 85 and its proof). 



C. Global and Eisenbud-Wigner time delay 

Thoughout this subsection we assume that, for some 9 > 4, tp is a wave packet in P™ 
such that Sip belongs to Vg nt . Then, as seen in the preceding subsection, one can study 
the limit of the local time delays T^(tp) and T^ ut (p) as R — > 00 by considering the limit of 
the quantities cf 1 r(p>) and a'^ lt (p), see (|^3*|) - (|^jl . Using the identity (St)* Si = U+ one can 
rewrite cr l ^(p) as follows (<p £ <pt = e~ lHlt p): 

POO 

= / (S m \[xt- R , R) ,St]<Pt) dt . (98) 
Jo 

Next, writing explicitly St = Stt + S r t and using the orthogonality relations (foTj) as well as 
the intertwining relations e~ lHat S a b = S^e - 6< one obtains 

POO 

<?%&) = I (StM[e iHot X(-R,R)e- iHo \Stt)tp)dt 
00 

iHot —iHot 







+ / (SrMie^x^n^e-^^S^dt 
Jo 

poo 

+ / (e~ iH ^Sttp\x(-R,R)e- lHrt S re p)dt 
Jo 

POO 

+ / (e- iH * t S r i<p\X{-Rjt>e- iH < t Su<p)dt. (99) 
Jo 

In the limit R — > 00 the last two integrals in ()99j) vanish. Indeed, since Sup £ 7Y_ and 
S r tp £ 7i+, the scalar products in the integrands converge to zero for each t £ R as R — > 00, 
and it suffices to justify the interchange of the limit and the integration. For this it is enough 
to majorize the absolute value of the integrands by an ^-independent integrable function. 
One has 



\{e-^Sttp\x W )e- tHrt S r tp)\< 



00 \ 1/2 / 1*0 \ 1/2' 

\(e- iH ^Sttp)(x)\ 2 dxj + U \(e- iH ^S rt p>)(x)\ 2 dx 



.(100) 
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By (jT§j) and (fT7|) with x = and 9 > 2, the r.h.s. of (jlUOj) is an integrable function of t on 
[0,oo). 

To treat the first two integrals in (|9*H|) we use the following asymptotic expression from 
Q: 

e^V^^' dt = f # ~ 1/2 " ^ + o(l) (101) 

as i? — > oo. This relation holds on D™ an d on T>2 Ut . The derivative d/dE in the second term 
on the r.h.s. means differentiation with respect to the energy variable E in Tt m as well as in 
H out (defined by flSEj-ffl}) H 

We first consider the terms proportional to that are obtained when applying to 

_____ _ 1 ^ 

the first two integrals in . Since S"« commutes with i/ , there is no contribution from 
the first of these integrals. By using the intertwining relation H r S r i = S r iHi one obtains 
the following contribution from the second integral in (|99|): 

^(S ri cp\{Ho 1/2 ~ (H + V r - V e r l / 2 }S re( p) . (102) 

IfVt^Vr and S r e<p ^ 0, the scalar product in ()102j) is strictly positive. Thus, ifVi<V r and 
tp is an incoming state (from the left) which is not entirely reflected by the potential, the 
local time delay t]£((/?) will not admit a finite limit as R — > oo, more precisely T^(ip) — > +oo 
as R — > oo. 

Next we observe that the terms from the first two integrals in (j99jl that are independent 
of R are just 

(Sutp\[-i-£Q,Stt]<p) + (S r e(p\[-iJg,S r e]<p) 



oo 



-dS u {E) 



iME)^^\^(E)dE 



POO 

/ ^(E)T ie (E)ip hl (E)dE , (103) 

JVt 

where %e(E) (the expression in the curly bracket) is one of the elements of the Eisenbud- 
Wigner time delay matrix T(E) at energy E which will be discussed below. 

A similar analysis can be carried through for <r^ ut (<£>), defined in (|96|) . using the asymptotic 
expression 

e iHot X { - R , R) e- iHot dt = -^H- 1 ' 2 - i^L + (i) . (104) 
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as R — > oo. One finds that 

D roo 

= --(S r M[Ho 1/2 -(H + V r -V e r 1 / 2 ]S reV )+ / ^(E)T u (E) V in (E)dE + o(l) . 
1 Jv t 

(105) 

Again, ii Ve ^ V r and if 99 is not completely reflected, T R nt (ip) will diverge as R — > 00, viz., 
r^ ut (</?) — > —00 as i? — > 00. However, the divergent term in a R nt (ip) is identical with that in 
cr R (<p), except for its sign. Hence the average <r R (<p) = | [cr^(^) + a R nt ((p)~\ converges to a 
finite limit, given by the r.h.s. of (|103|) and denoted by r((p). 

The preceding result can be rewritten in terms of the Eisenbud-Wigner time delay oper- 
ator T. In the representation 7i m of the Hilbert space 7i (see f!66j) ) T acts at energy E as 
an operator T(E) given as follows: 



(i) for V t <E< V r , T(E) is just multiplication by T U (E) = -iS u {E) dS ee (E)/dE, 

[ii) for E > V r , T(E) acts on tp m (E) 2 matrix, i.e., 

T(£) = 1^ Tfr(E) I , (106) 



T H {E) T rr (E) 



where (for a, b = i or r) 



T ab (£) = -i £ SUE) = -i £ S*(E) ac ^l. (107) 

c=£,r c=£,r 

Since the incoming wave functions considered in this subsection have only one non-zero 
component (v?J. n (i?) = 0), the r.h.s. of ()1U3|) is just the mean value of the operator T for the 
initial wave packet tp: 

r(<p) = lim tM = lim 1 [t£(^) + r° ut (tp)] = (<p\T<p) . (108) 

If /i > 2, T(_E) is well defined for E ^ Vg, V r (see Appendix D). Also, as a consequence 
of the unitarity of S(E), the matrix T(E) in ()106|) is hermitian, and T^(E) is real for each 
E > Ve (E ^ V r ). This implies that the family {T(E)} determines a unique (in general 
unbounded) self-adjoint operator which we have denoted by T. This time delay operator T 
commutes with H m and constitutes a quantum-mechanical observable. The next subsection 
is devoted to some further comments on the meaning of this operator. 

To end this subsection we point out an interesting alternative expression for the global 
time delay r(<p). We observe that, for <p G TC + as well as for (j) G H- and any k G R, one 
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has the following identity: 



00 J Xl 



By using this identity in T R n (ip) and T^ ut ((/?) one may write T R ((p) as 22| 



(109) 



r R (ip) = T R/ (ip) + T R>r (ip) + as R -> 00 , (110) 



where 



r>00 /"0 

= / dt / dx {|(e- iH ^7^(x)| 2 -|(e-^V)(^)| 2 -|(e"^^)(^)| 2 } 



00 J-R 
00 /-ij 



t r M = / dt / dx {\(e- im njif)(x)\ 2 -\(e- iH ^S r ^)(x)\ 2 } . 

J-oo Jo 

By arguments similar to those applied before one finds that both r R) e((p) and r R ^ r (ip) 
converge as R — > 00, their limits (denoted Tg(<p) and r r ((p) respectively) being given by 
setting R = 00 in the double integrals defining T R> e(<p) and r R>r {ip). Thus 

rfa) = r t (ip) + r r (<^) . (Ill) 

This shows that the (symmetrized) global time delay r((p) may be decomposed into a con- 
tribution associated to the left spatial half-interval (—00, 0) and a contribution coming from 
the right half-line (0, 00). One observes that the final state Sg(p appears in Ti((p) only through 
its reflected part S^ip and in r r (ip) only through its transmitted part S r i<p. 



D. Discussion 



In the literature and in text books one usually uses the Eisenbud-Wigner expression 

(<p\T<p) = ^^{-z^^^-i^^5}^(E)dE (112) 

to compute the time delay induced by a one-dimensional potential V on an incoming wave 
packet tp G 7i + . Its identity with the symmetrized global time delay, Eq. (|108[) . allows one 
to have a time- dependent interpretation of what is computed with (|112jh An alternative 
interpretation is furnished by the decomposition (jlllj) of the symmetrized global time delay. 
The symmetrized local time delay T R ((p) can be written in the following form: 

tr&) = \[tM + < (*>)] = T R {nj<p) ~ I [T'M + IT (*>)] ■ (H3) 
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We have shown that tr(<^) is the appropriate expression admitting a finite limit as R — > oo. 
This implies that the pertinent reference sojourn time T^ f ((p) is neither T^((p) nor T^ ut (y?) 
but their average, i.e., 

T^) = \[T^)+Tr^)] ■ (H4) 

Notice that T^ { (ip) depends on the potential V since the free outgoing sojourn time 
Tft nt ((p) involves the final state Sip. There are however two special cases where T^ { (ip) is 
actually independent of V: 

(i) The first one is the situation where Vg = V r . In this case the difference T l ^{ip) — T^ xt {ip) 
vanishes as R — > oo, so that T]f(ip) = T^(ip) + o(l) as R —> oo. 

(ii) The second one is the general situation < V r but with an incoming wave packet 
ip having energy support (relative to Hi) contained in the interval (Vg,V r ). Since such an 
incoming state is completely reflected by the scatterer V one has \(Suip)(—p)\ = \0(p)\ for 
all p, so that T% ut (ip) = T%{ip) and therefore T T R ei (ip) = Tg 1 (</?). 

The statement in (i) is easily obtained by using (jlOlj) and (|l()4|l . that in (ii) is a consequence 
of pity . 

In the above two cases any of the incoming, outgoing and average reference sojourn time 
is acceptable, i.e., one has in these cases: 

r(ip) = lim T%{ip) = lim r°^(ip) = lim r R { V ) . 

H. — >oo H — >oo H. — >oo 

Therefore one can interpret the Eisenbud-Wigner expression (jll2}) as the difference between 
the sojourn time Tn(Qjip) (where Qj(p evolves with H = H + V) and the free incoming 
sojourn time T^(ip) (where ip evolves with Hp = Hq + Ye) as R — >• oo. 

We next discuss the situation where Vi < V r and tp G TL+ is an incoming wave packet 
having energy support (relative to He) above the threshold V r (i.e., in (V r ,oo)). Here the 
reference sojourn time ((f) will depend on V (and not just on its asymptotic values Ve 
and V r ), so that one cannot have a similar explanation of (jll2J) as the one obtained for the 
cases (i)-(ii) above. In this case one can still use any of the two expressions given in (jll3|) or 
the alternative expression (jlllj) as a time-dependent interpretation of ((p\T<p). Nevertheless 
it is instructive to introduce the following reference potential: V = VeXe + VrXr (i-e., the 
step-potential represented in dotted lines in Fig. 1). Notice that V depends on V only 
through its asymptotic limits Vi and V r , and that the point of discontinuity of V coincides 
with the center of the interval [-R, R] used to define the local time delay Tr((p). 
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For the step-potential V one knows explicit expressions for the quantities of interest in 
scattering theory (see e.g., Chapter 1 in |35j]). For E > V r all entries of the S-matrix S(E) 
are real. It then follows from the unitarity relation for S(E) that the diagonal elements 
(but not the off-diagonal ones) of T{E) are zero. Therefore the Eisenbud-Wigner expression 
(I112j) vanishes if (p G Ti. + has energy support (with respect to Hg) above V r In this 

situation ()112|) is non-zero only if the potential V is different from V. In other words, 
only the "wavy" part V — V of the potential V may induce a time delay on tp and one 
may therefore interpret the Eisenbud-Wigner expression as the effect of V — V. Note that, 
although the scatterer V induces no time delay above V r , it still affects the incoming wave 
packet <p (S I ii H = Hq + V) . Note also that any potential V with Vg < V r induces an 
infinite incoming global time delay above V r (r in (<^) = +00) so that the Eisenbud-Wigner 
formula (or equivalently the symmetrized global time delay) is the appropriate expression 
to measure the finite effect of the wavy part V — V of the potential V. 

In conclusion, for a state if incident from the left, one may interpret the symmetrized 
global time delay, or equivalently the Eisenbud-Wigner expression, as the effect of the full 
potential V on the components ip m (E) of (p with energy E below V r and as the effect of the 
wavy part V — V of V on the components ip m (E) with energy above V r . 



VI. CONCLUDING REMARKS 



We have presented a general scattering theory and obtained the existence of the global 
time delay as well as its identity with the Eisenbud-Wigner expression for potentials V 
having different limits at x — —00 and at x — +00, assuming V to be bounded and 
approaching its limits at x = ±00 at a certain minimal rate (specified by the number ^ in 
(P"U|) - fTT|) ). From this identity we obtained a time- dependent interpretation of the Eisenbud- 
Wigner expression. The same results can be established under weaker assumptions on V by 
using more refined versions of our approach or different techniques. It is possible to handle 
potentials with (square-integrable) local singularities and to weaken the assumptions that 
we made on /i. The hypothesis /i > 1 (short range condition) is sufficient for obtaining all 
results on scattering theory in Section HVl j^. As regards time delay we expect (in view 
of known results |3i| for n-dimensional Hamiltonians with potentials converging to zero 
as \x\ — > 00) that the symmetrized global time delay should exist and coincide with the 
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Eisenbud-Wigner expression if fi > 2 and for initial states belonging to T>™ for some 9 > 2. 

Only wave functions with energy support away from the thresholds Vi and V r have been 
considered. At thresholds the time delay is usually infinite. As a consequence the Eisenbud- 
Wigner operator T will be unbounded, and its behavior near the thresholds would require 
a more refined investigation. 

In one-dimensional scattering problems it is interesting to distinguish between the re- 
flected and the transmitted part of a wave function, and there is a considerable literature 
on tunnehng times <Q, and references cited therein). The va.ne of the g .o b a. ti,ne 
delay for a given initial state involves both the reflected and the transmitted wave and is 
therefore not related directly to a tunneling time. However the global time delay has the 
merit of being given in terms of a self-adjoint linear operator and can thus be interpreted 
as a quantum-mechanical observable in the usual sense; furthermore it has a meaning also 
for scattering systems in more than one space dimension. 

We point out that no oscillating terms in R were involved when we considered the limit 
of the local time delays as R — > oo. In various other publications on time delay (e.g., Iljl, 

n n n 

|15| . |41|, t ne authors encountered oscillatory terms like sm(2pR) and then presented an 
argument to suggest that these terms will not contribute to the global time delay. In it 
is stated that these oscillatory terms are related to the uncertainty principle. Actually, the 
presence of such terms arises when one works with non-normalizable eigenfunctions of the 
Hamiltonian. In a fully Hilbert space derivation of global time delay, with sqnare-integrable 
wave functions, one has no problem with oscillatory terms (see also |4l|). 

One could also consider the global time delay r Xo (<p) obtained by starting with the local 
time delay T^°((p) in the translated interval [— R + Xq, R + xo], where Xq G M. It is clear that, 
for tp G H + , the local time delay t^°(<£>) is given by (0) with the following substitutions: 



(f i — > ip x ° = e lPx °(p and S i— > S x ° = e lPx °Se lPx °. Then, proceeding as in 42 1 or |43[, one 
obtains 

r*%) = rfa) - ^{S u <p\[p-\S u ]y) - ^{S t M[P-\Sh]<p) ■ (115) 

Thus two additional terms appear in the above situation. The presence of such terms was 
already pointed out in the case Vg = V r j3| and also in a more general situation [42J]. 
Recalling the identity (jl()8j) one sees from that the Eisenbud-Wigner expression for 

time delay assumes implicitly that the spatial interval [-R + Xo, R + Xq] on which the total 
and reference sojourn times are compared is centered at the origin (xq = 0). Note in 
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particular that the relation (|115|) implies that the time delay r(ip) due to the translated 
step-potential V = V^X(-oo,x ) + KX(a; ,cc) (^o 7^ 0) for a wave packet </? G 7i + having 
energy support (with respect to He) above V r is non-zero in general. Finally, notice that the 
sojourn times ()83j) - ()85)) and therefore the symmetrized time delay r((p) are invariant under 
time translations. 

As a last point we mention that one can find a general study and discussion of the 
symmetrized definition of time delay in (3). These authors point out in particular that the 
symmetrized time delay is useful in multichannel scattering processes, invariant under an 
appropriate mapping of time reversal and relatively insensitive to the shape of the spatial 
regions used for defining the local time delay in more than one space dimension. 

ACKNOWLEDGMENTS 

This investigation was suggested to us by S. Richard. We also have benefited from 
numerous discussions with him and with V. Cachia and R. Tiedra de Aldecoa to all of 
whom we wish to express our sincere gratitude. This work was partially supported by the 
Swiss National Science Foundation. 

APPENDIX A: TIME DECAY IN A CONSTANT POTENTIAL 

We show here how to obtain (fTTj) and (|2Tj). The estimates (fTHjl - PUj) can be deduced in 
a similar manner. We assume (without loss of generality) that xo = in ()17)) . We first 
consider the case \t\ < 1 and let tp be any wave packet. If g is a bounded function, say 
\g(x)\ < M < oo, then for \t\ < 1 and 9 > 0: 



Taking g satisfying g{x) = 1 if x < and g(x) = otherwise, or g = f with / as in (|21j) . 
one sees that the inequalities (j!7)) and (|?H) are satisfied for |t| < 1. 

To treat the remaining values of t (t > 1 in ()17|) and \t\ > 1 in (j21j0 . we write 9 in the 
form 9 = In + 2e, with n a non-negative integer and < e < 1, and we shall use Taylor's 
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formula for the function F(s) = e ts (s G M), i.e., 

k=0 

with 



|i«n(s)| < (1 + e)|sT +e <A\s\ n+e (A-2) 

for < e < 1. 

Let ip be a wave packet with positive momentum and let g = 1—h, where h is the Heaviside 
function (so g(x) = for x > and = 1 for x < 0). Let t > 1. Then g(2tp) = for 
all p > 0, so that g(2tp)0(p) = 0. Also, since Q is just differentiation in momentum space, 
the wave function Q k ip has the same momentum support as ip for any positive integer k 
(assuming at least k times differentiable). Hence g(2tp)(J-'Q k (p)(p) = i k g(2tp)0^ k \p) = 
for all real p. Thus, taking for example n — 1 (i.e., 2 < 9 < 4), we get from (fTHj) that 

poo 

\<p t (x)\ 2 dx = / | fl (2tp)[^/\)](p)| J dp 



00 



00 q2 

\g{2tp)[T{e^ 2 '\ -<p- &-<p)\{p)\ 2 dp • 
■00 41 

By using the bound |p(2£p)| < 1, the unitarity of T and (|A-1|) - (|A-2|) with n — 1, one obtains 
that 

|^(x)| 2 rfa; < / \(e ix ^-l-t-Mx)\ 2 dx 
■00 J —00 ™ 

1 f 00 

< Jr^)/ IM 2(1+e V(*)l 2 ^ 

which implies the validity of (fTTjl for t > 1 and # G [2, 4). The result for other values of 9 is 
obtained similarly (using an n-th order Taylor expansion if 9 = 2(n + e)). 

We finally show how to verify ()21|) for \t\ > 1 and 9 e [4,6). Let / be as stated in (c). 
Using the inequality \a + (3\ 2 < 2\a\ 2 + 2\/3\ 2 we get from ([16]) that 



|/(x)^(x)| 2 rfx = / |/(2tp)| 2 |[^(e^/ 4 V)](p)| 2 ^ 

/°° O 2 O 4 

|/(2tp)| 2 |[^ 2 / 4 V - ¥> - ^ + ^^)](P)| 2 

+ 2 y |/(JMp)| a |[^ + i^-^)](p)|>. (A-3) 
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In the first integral on the r.h.s. of (|A-3|) we majorize \f(2tp)\ 2 by C 2 and use Taylor's 
formula (with n = 2) to obtain, as above, the following upper bound for this term: 



which is majorized by the first contribution in (|2~Tf . with 9 = 2(2 + e) (and Cg = 2 e+1 C 2 ). 
To treat the second term on the r.h.s. of ()A-3|) we set rj — </? + iQ 2 (p/4:t — Q A (p/32t 2 and 
observe that fj(p) ^ only for \p\ > p if if has momentum in R \ (— Po,Po). We majorize 
\f(2tp)\ by C(l + 2|p ||t|)~ M in the integral and get 



which is majorized by the second contribution in ()21)1 if \t\ > 1. 
APPENDIX B: THE MOURRE ESTIMATE 

We establish here the validity of a Mourre estimate, as stated in Section III Bl for our 
class of Hamiltonians (/i > 1). We shall freely use the following properties of compact 
operators. The product of a compact operator and a bounded operator is compact. If {K n } 
is a sequence of compact operators that converges in operator norm, i.e., if there exists a 
bounded operator K such that \\K n — K\\q^ — > as n — > oo, then K is also compact. If f 
is a bounded function on 1R satisfying f(x) = near x = ±oo or, more generally, such that 
f(x) — > as \x\ — ► oo then, for our class of Hamiltonians, the operator (H — z)~ Y f{Q) is 
compact for each non-real z. Also, if F(A) denotes the projection onto the subspace of wave 
functions having energy support (with respect to H) in the interval A, then f(Q)F(A) and 
f(Q)PF(A) are compact if / is as above and A is a bounded interval. 

We fix a smooth function J r satisfying < J r (x) < 1 for all x, J r (x) = for x < 
and J r (x) = 1 for x > 1, and we set Jg = 1 — J r . We denote also by J r the operator of 
multiplication by J r (x). The following decomposition of V will be used: 



We consider wave packets i[) with energy support in a finite interval A = [a, /3], i.e., satisfying 
vj) = F(A)ip. For A = (PQ + QP)/A, one has [iH ,A] = [iP 2 ,A] = P 2 = H - V and 
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V — (V - V e )J e +{V- V r )J r + V t Ji + V r J r . 



(B-l) 



[iV, A] = —QV/2 (assuming that V is differentiable, see the Remark below). Hence 
MiH,A]^) = ^m)-{^\V^)-^\QV'i>) 

-(^\(V - Vt)W) - mv - V r )J r ip) - ^iQVty) • (B-2) 

(i) We first assume that V r < a < j3 < oo. One has 

v t (ii>\Jd>) + v r {if>\j r ij>) < VM{J* + <W> = vm^) ■ 

Hence flB-2jl leads to 

A]ip) >(a- K)(V#> - (MV - V t )J t + (V- V r )J r + l -QV'}^) . 

Since (V(x) — Ve)j£(x) + (V(x) — V r )J r (x) + ^xV'(x) converges to zero as x — ► ±00 (provided 
that V' decays more rapidly than | a; | 1 ) , the Mourre inequality (}28j) holds with A = a — V r 
and K = -F(A)[(V - V e )J e + (V - V r )J r + \QV']F{A). 

Remark. If V is not differentiable the compactness of F(A)[iV, A]F(A) is obtained (for 
fi > 1) by decomposing V as in (jB-l|) and writing for example 2[iViJi,A] = —ViQJ[ and 
2[i(V - V e )J e , A] = i[Q(V - V e )J e ]P - iP[Q(V - V e )J e ] + (V - V t )J e . 

(ii) We now consider the case where < a < (3 < V r . Since we assumed that ip = F(A)ip, 
the sum of the last three terms in ()B-2j) is again of the form {i/j\Kqi(j) for some compact 
operator Kq. We shall show that, for A C (Ve, V r ), the operator J r F(A) is compact. Writing 
Ve(ij)\Jeip) = Vg(ip\ijj) — Vi{ijj\J r ip), it then follows from ()B-2j) that a Mourre estimate holds 
with X — a — Vi and K — K + F(A)(V t - V r )J r F(A). 

Let g be a smooth function defined on K such that g(E) = 1 for E G A, g(E) = for 
E > (/3 + V r )/2 and for < Vg. If ^ has energy support in A, then ^ = g(H)tp, which 
implies that g(H)F(A) = F(A), and it suffices to show that J r g(H) is a compact operator. 
For this we introduce an auxiliary Hamiltonian H = Hq + V, with V(x) = V(x) for x > 
and V^x) = V r for x < 0. We may write 

J r g(H) = J r g(H) + J r [g{H) - g(H)] . (B-3) 

The operator g(H) (hence also J r g(H)) is compact (even of finite rank), because the Hamil- 
tonian H has only discrete spectrum below its threshold V r , hence at most a finite number 
of eigenvalues below {(3 + V r )/2. 

40 



To handle the second term on the r.h.s. of (|B-3|) . we use the following formula for g(H) 



(Theorem 6.1.4 and Remark 6.1.3 of 271]): 



- (H-E + i)- 1 } dE 



9{H) = in £ 9(E)[(H-E-i 

+ — / g'(E) [(H-E- I)" 1 + (H-E + s)" 1 ] dE 

1 /»1 POO 

-— ede g"(E)[(H - E - ie)' 1 - (H - E + ie)~ l }dE . (B-4) 
2tt2 Jo y_oo 



The integrals in (|B-4|) exist in operator norm (i.e., the approximating Riemann sums converge 
in operator norm). By using (|B-4|) and the corresponding formula for g(H), one has (with 
the notations R z = (H — z)~ l and R z = (H — z)' 1 ): 

1 f°° 

J r [g(H) - g{H)} = — / g(E) [J r (R E +i - R E+ i) - J r (R E -i - R E -i)\ dE 

2tt? y.oo 

+ — g\E) [J r (R E+i — R E+ i) + J r (R E -i — R E -i)} dE 

rl poo 

— ede g"(E)[J r (R E+ie - R E+l£ ) - J r (R E - ie - R E _ le )] dE . (B-5) 

=T« JO J-oo 



2m 

Since the integrals in ()B-5|) exist in norm, this implies the compactness of J r [g(H) — g( H)] by 
using the fact (established below) that J r (R z — R z ) is a compact operator for each non-real z. 
To verify the compactness of J r {R z — R z ), we write 

J r (R z -R z ) = J r [(H - z)- 1 - (H - z)- 1 ] 

= J r (H - z)- l {V -V)(H - z)~ l 
= (H - z)- l J r {V -V){H - z)- 1 + 

+(H - Z y l [H - z, J r }(H - zy\V - V)(H - zY 1 . (B-6) 

The first term on the r.h.s. of (|B-6)1 is zero since J r (x)[V(x) — V(x)} = 0. In the second term 
we observe that [H — z, J r ] = [P 2 , J r ] = — J" — 2iJ' r P, and the compactness of the second 
term follows since J"{x) = J' r {x) = for \x\ > 1 and P(H — z)" 1 is a bounded operator. 

APPENDIX C: EXISTENCE OF LARGE TIME LIMITS 

In this appendix we prove the existence of the (strong) limits involved in the definitions 
of the scattering projections Ff and Ff- (see (|4T |) - (jl2*|l ) and of the operator W given by 
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Eq. f)6Uj) . It suffices to establish the existence of these limits on a dense set of wave functions 
in TC C (H). We consider the dense set of if} having bounded energy support (with respect to 
H) away from the thresholds Ve and V r . 

(a) We first consider F^~, as given by (|4fi|) . for which we use (f3*H|) with h,\ = hi = H, hence 
W T = e tHT g(Q)e~ lHT . Due to the first projection F C (H) in (|4fij). the supremum over the set 
{(f> G H | ||0|| = 1} in can be replaced by that over the set {0 G TC C (H) \ \\<f>\\ = 1}. By 
the Cauchy-Schwarz inequality the integrand in (J33J), with Vi = V2, is majorized by 

||(1 + IQD^e-^H • || [(1 + \Q\)g"{Q) + 2i(l + |Q|)</(Q)P] e -^|| . (C-l) 

Inserting this bound into (|35j) and then applying the Schwarz inequality in L 2 ([s,t}) (i.e., to 
the integral over the interval [s, i]), one has 

pt poo pt 

\\Wti)-W 8 ip\\ 2 < sup du dx |[(l + |Q|)-V^](x)| 2 - / N 2 dr, (C-2) 

4>&H C {H), ||0||=1 is i-00 is 

where A^ T is defined as the second factor in (jC-lj) . Let E be a bounded closed set in (Ve, 00) 
disjoint from V r . By (13 lj) there is a constant Cs such that 

/oo 
dxita+iQirv^Kz)! 2 ^^ 11011 2 
-00 

for all <p in 7i c (H) having energy support in E. By restricting the supremum in ()C-2Jl to 
this set of 0, one finds that 

\\F(E)W t ip - F(E)W s iP\\ 2 <C E f N 2 dr, (C-3) 

J s 

where F(H) denotes the projection onto the subspace of wave functions having energy sup- 
port (with respect to H) in E. We show below that J °° N 2 dr < 00. Hence the sequence 
{F(Y?)W T ip} is Cauchy as r —> +00 and therefore s-\im t _ >+OQ F(E)e lHt g(Q)e~ lHt F c (H) exists. 
By varying the set E, one can conclude that s-lim t ^ +00 F c (H)e lHt g(Q)e~ lHt F c (H) exists Q|. 

We now comment on the finiteness of J °° N 2 dr. For the first term in N T (i.e., ||(1 + 
\Q\)g"{Q)e~ lHT ip\\) this is immediate from (|3*Tj) since g"(x) = outside the interval [—1,0], 
hence (1 + \x\)g"(x) < C(l + |a;|) _1 . The second term in iV r requires more care because P is 
an unbounded operator and does not commute with e~ lHt . By a commutator calculation (see 
(c) below) this term can be expressed as a sum of four terms of the form Df(Q) (H + i) m ijj T , 
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with D a bounded operator, / satisfying \f(x)\ < (7(1 + \x\) 1 and m — or m — 1. Clearly 

/oo /»oo 
rfr / dx \[Df(Q)(H + i)^ T ](x)\ 2 
-oo J — oo 

/oo 
11(1 + |Q|)-V^(tf + zy>]fdr < oo (CM) 
-oo 

since (H + belongs to 7i c (H) and has the same (bounded) energy support as ip. 

(b) The arguments for the existence of s-lim t _ >+OQ e lHet g(Q)e" lHt F C (H) (implying that of 
W in (|60[)) are very similar to those used in (a) above. The only difference arises through 
the fact that now hi = (instead of hi = H). This leads to an additional term in J N 2 dr, 
viz., 

pt POO 

dr dx |[(1 + \Q\)(V - V t )g{Q)e- iH ^]{x)\ 2 . (C-5) 

J s J — oo 

Again it suffices to know that the integral with respect to dr in (|C-5|) is finite if t = +00; 
this follows from flHU) since |(1 + \x\)[V{x) - V e }g(x)\ < C(l + l^l)" 1 if // > 2 in (fTU|) . 

(c) In (a) above we used the following formula, with <p = (1 + I^Dfi 1 ' a smooth function 
vanishing outside some finite interval A: 

4 

0(Q)P = ]T D k f k (Q)(H + *H (C-6) 
fc=i 

with Dk bounded operators, f k smooth functions vanishing outside A and m k = or 1. This 
can be obtain as follows: 

<P(Q)P = 4>{Q)P{H + i)-\H + i) 

= P(H + z)~V(Q)(# + + [<P(Q), P{H + *) _1 ](# + z) . 
Evaluation of the commutator gives 

[0(g), p(H + iy l ] = [0(g),p](F + z)- 1 + p[0(g),(i7 + z)- 1 ] 

= i(f)'(Q)(H + zy 1 - P(H + zy'l^Q), H + i](H + i)" 1 

= i<j)'(Q)(H + z)" 1 - P(H + z)~ 1 {-0"(g) + 2iP<P'(Q)}(H + z)" 1 . 

So 

<P(Q)P = P{H + i)- l ^{Q){H + i)+i(j}'{Q) 

+P{H + i)~ 1 (j>"{Q) - 2iP(H + i)- 1 P<j ) , {Q) , 

which is of the form (jU-6|) since P(H + and P(iJ + z) _1 P are bounded operators. 
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APPENDIX D: DIFFERENTIABILITY OF THE S-MATRIX 



The entries of the 5"-matrix S(E), i.e., the transmission and reflection amplitudes at 
energy E, are given as simple functions of E combined with Wronskians of Jost solutions for 
the potential V. We refer to or j3] for a complete description of these expressions and 
to Section XVII. 1 of Q for a presentation of their derivation. For example the transmission 
amplitude S r g(E) at energy E > V r is just 

S " {E) - wumim • (D_1) 

where kg. — (E — Vj?) 1 / 2 , k r — (E — Vr) 1 / 2 and the Wronskian of the left and right Jost 
solutions fi and f r is (' = d/dx) 

W{h{E)J r {E)) = f t (E,x)ft(E,x) - f r (E,x)ti(E,x) . 

Thus differentiability or Holder continuity of the Jost solutions and their spatial derivatives, 
as functions of E, imply differentiability or Holder continuity of S r g(E). 

Properties of this type can be obtained from the Volterra integral equations for f% and 
f r . Consider the following integral equation for a function ^f KjCr (x), where k and a are two 
parameters: 

V^ix) = $(«, a- x) + ^ SmK{y ~ X) V(y)y K Ay) dy , (D-2) 

Jx K 

assuming that the inhomogeneous term ct; x) satisfies an inequality of the form 

|$(k, ct; x)\ < ti(k, a)h(x) (D-3) 



with h(x) > 1 and 

POO 

h(y)\V(y)\dy < oo . (D-4) 



o 



The standard iterative method for solving Volterra equations gives the following bound on 
the solution ty K;(T (x) for x > 0: 

\^ K A^\<V^^)Kx)exp(^J^ h(y)\V(y)\dy^ . (D-5) 

The integral equation for the Jost solution f r is of the above form, with \l/ Kj(T (x) = *$? K (x) = 
f r (E,x), k = k r = (E - V r ) 1/2 , V(x) = V(x) - V r and $(k,ct;x) = e iKX , hence t)(k,(t) = 1 
and h = 1. The existence of f r (for E > V r ) requires the assumption that fi > 1 in (fTTjl . The 
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derivatives /i"' of f r with respect to k also satisfy an integral equation of the form (jD-2|) 
obtained by formally differentiating the integral equation for f r and calculating at each step 
the inhomogeneous term $(k,<j;x) = <&(k,x) (the parameter a plays no role here). By 



considering successively / r , fr 2 \ . . . and using at each step the bounds of the form (|D-5|) 
obtained for the derivatives of the lower orders, one finds that the inhomogeneous term in 
the integral equation for fr satisfies (|D-3|) with ij(k, <j) = c(k) and h{x) = 1 + x n (x > 0). 
The number c(ac) is finite if k > and y n \V(y) — V r \ dy < oo. So the integral equation 
for f^ has a unique solution provided that /i > n + 1 in By using the bounds for the 

functions / r (0 < j < n), one can also verify a posteriori, proceeding again recursively, that 
these functions are indeed the derivatives of the Jost solution (using a theorem permitting 
the interchange of the integral with the derivatives with respect to k in the occuring Volterra 
equations, e.g., Lemma 2 in Chapter XIV of 4^|). 

Next let k, a > with \k — a\ < 1 and let \l/ Kj(T be the difference of the Jost solutions 
f r at energy E K and E a : \1/ K>0 - = f r {E K ) — f r (E a ), where for example E K = k 2 + V r . The 



function ^ KiC satisfies an integral equation of the type (|D-2J) . By using the inequality 
\ e ^-e iax \ < 2\k — <j\ i/ \x\' / valid for < v < 1, the addition theorem for the sine function and 
the bound ()D-5j) for f r (E a ), one obtains an estimate of the form ()D-3)1 for the inhomogeneous 
term $(k, a; x), with r)(K, a) = p(n, cr)\n — cr\ v and h(x) = 1 + \x\ u , where p(n, a) is a smooth 
function (the condition (jD-4|) is satisfied if v < fi — 1). In view of (|D-5|) we conclude 
that f r (E) is Holder continuous as a function of E > V r with any exponent v satisfying 
v < min{l,/i — 1}. By proceeding again recursively, one can use the integral equation 
for fr H \E K ) — fr n \E a ) to show that fr H \E) is Holder continuous with any exponent 7 < 
min{l, ji — n — l}if//>n+l. 
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